The superheated Melting of Grain Boundary 
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Based on a model of the melting of Grain Boundary (GB), we discuss the possibility of the 
existence of superheated GB state. A Molecular Dynamics simulation presented here shows that 
the superheated GB state can realized in the high symmetric tilt GB. Whether the sizes of liquid 
nuclei exceed a critical size determined the superheating grain boundary melting or not. Our results 
also indicate that the increase of melting point due to pressure is smaller than the superheating due 
to nucleation mechanism. 

PACS numbers: 68.35.Rh,64.70.Dv,68.35.Fx,68.35.Ct 
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I. INTRODUCTION 

The superheating has been found in a larger number 
of systems such as surface QHQ , small cluster [4|, 
confined thin film [f| and particles covered (or embedded 
in) by material with higher melting point 0. Gener- 
ally, the melting of solid material is heterogeneous pro- 
cess with the nucleation mechanism at surfaces or inter- 
faces 0, H| ■ Providing heterogeneous nucleation could be 
avoided by means of suitable coating 6] or internal heat- 
ing Q , the metal crystal can be in the superheated state, 
its melting is completed by a thermodynamically insta- 
bility resulting in homogeneous disordering and catas- 
trophic mechanism with the stability limit from 0.2T m 

to2.or m d miii mm. 

A large number of researches have been contributed to 
the role of surface for the melting of crystal. The super- 
heated melting of fc c(llO) a nd fcc(lOO) surfaces are virtu- 
ally never observed 0, 0, 0] . The only example of the 
superheated surface is the small crystal strictly confined 
by high-symmetry fcc(lll) facts [HHl-III]- GB as another 
important quasi-2D defect also leads to the heteroge- 
neous melting of solid material. Quite a number of stud- 
ies have shown that GB can't melt below the T m (not pre- 
melting) [SGSIIIIIJIIIIIIIMI^ Us- 
ing the MD simulation, Kikuchi and Cahn, Ciccotti et.al. 
showed that GB doesn't melt until temperature reaches 
to melting point of bulk. Nguyen et.al. using the more 
accurate interatomic potential by embedded-atom meth- 
ods(EAM) studied the high-temperature GB structure, 
they found that, close to melting point T m , the GB struc- 
ture was disordered, quite liquid-like and meta-stable, 
and over a long interval of simulation the underlying crys- 
talline order can re-emerge. The experiment (T.E.Hsich 
and R. B. Balluffi) using the HREM(High Resolution 
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Electron Microscopy) methods supported above argu- 
ments and showed that aluminum GB did not melt below 

o.999T m [ail- 

It is more surprising that some the simulations also 
hint that some high symmetric GBs similar to high sym- 
metric surface can probably melt by the superheated |24j . 
In this work, parallel to Di Tolla's work for surface 
we study the possibility of superheated high-symmetry 
GB by a theoretical model and MD simulation of a sym- 
metric aluminum GB. 

When temperature beyond melting point a crystal 
reaches the superheated state, all liquid nuclei must be 
smaller than a critical size. These liquid nuclei are un- 
stable and able to re-crystallize again. The liquid nuclei 
easily form at surfaces, grain boundary and other solid 
defects regions. So, to avoid larger liquid nucleus, the 
crystal must be prepared with the lowest number of solid 
defects. It's easier to study the superheating of crystal 
in computer simulation than in experiment. We can con- 
struct prefect crystal using the periodic boundary condi- 
tion in computer simulation. In experiment it's very dif- 
ficult to obtain infinite volume prefect crystal. We can't 
eliminate the influence of surfaces, grain boundaries, dis- 
locations and other defects with complicated structures. 
All these defects have potentially become the liquid nuclei 
to melt crystal. However surface effects can partially re- 
move by coating other material with higher melting point 
or internally heating the material. By these methods we 
can obtain superheated crystalline grains. 

Additionally, the grain boundary itself may probably 
become the liquid nucleus near the melting point of crys- 
tal. Based on our theoretical model, it's possible that, 
under condition of the absent of critical nucleating cores, 
grain boundary doesn't melt even the temperature be- 
yond melting point. We hope to find the superheated 
grain boundary in computer simulation, although it is 
difficult to find in experiments due to different kind of 
unavoidable nucleation mechanism. Our simulation will 
show that high-symmetry grain boundary can sustain 
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the melting is generally completed by the mechanism of 
heterogeneous nucleation. We consider that a liquid film 
with thickness 21 forms between two semi-infinite solid 
(Fig. ma)). The change of free energy per unit area is 
taken as 



AF(l) = 2pLl(l - T/T m ) + A 7 (0 



(1) 



where p is the liquid density, L the latent heat of melt- 
ing, A 7 (Z) the difference between the overall free energy 
Isl-sl of two interacting solid-liquid interfaces sepa- 
rated by a distance I and the GB energy per unit area 



A j (I) = 7SL-SL - 7gb 



(2) 



By extending the Cahn's wetting theory to solid-solid 
interface [U, using A7(0) = and only considering the 
short range interaction, we may obtain 



A 7 (0=A 7oo (l-e- i /S). 



(3) 



where A 700 = 2^/sl—Igb is the difference of the interface 
energy of two isolated solid- liquid interface ^sl and the 
GB energy 7gb, £ is the width of solid-liquid interface. 

The condition of GB melting is defined by following 
process. We may image there is a droplet in GB region 
(Fig. [Hb)), the equilibrium condition of the droplet is 



FIG. 1: (a) The change of atomic density across the melt- 
ing GB. (b) The schematic of partly wetting GB. (c) Critical 
liquid thickness of a non-melting GB vs temperature above 
T m (schematic). Inset: free energy change upon conversion 
of a film of thickness I from solid to liquid. From T m to T 3 
the solid GB is stable. (d)The structure Al E=13 GB. The 
solid and open circles show the I and 11(001) planes of the ...I 
II I II..., stacking sequence along z. The basic GB structural 
units are also shown with dash-dot lines. This figure illus- 
trates only part of region near grain-boundary plane of our 
simulation cell 



above the melting point of crystal without other nucle- 
ation mechanism. The superheated grain boundary is 
easily understood by proximity effects if we consider the 
grain boundary is sandwiched between two superheated 
crystalline grains, while the superheated grains can be 
obtained by properly internal heating. 



II. THE MODEL OF GB MELTING 

The melting point T m of a solid may be defined as the 
temperature with the coexistence of solid phase and liq- 
uid phase. For a solid with surfaces or grain boundaries, 



7gb - 2jslcos(6) = (4) 

where 9 is the wetting angle. If A 700 > 0, i.e. job < 
27SL: & is a finite value, the droplet can survive in the 
region of GB (partial wetting), the whole GB can't be 
wetted. As Aj^ < 0, i.e. ,jgb > 2jsl, can't be 
defined. A liquid film forms in GB region (Wetting) and 
the GB melts. Thus A7oo = 0, i.e. jgb = 2^sl and 
6 = 0, may be considered as the criteria of GB Melting. 

In early time, one was decline to think that GBs melt 
below melting temperature, that is, at certain tempera- 
ture below T m , A7QO < 0. In this paper we only consider 
the possibility of GB superheated melting, that is, as 
T > T m , A loo > 0. 

For superheated melting, the insert of Fig. W[c) show 
that AF(l) has a local minimum at I — 0, and approaches 
to negative infinite as I approach infinite. There is a 
maximum at l c , which is 



l c = iLn{ "°° Jm - ) (5) 

The Fig^c) shows the temperature dependence of l c . 
According to Fig^c), at a certain temperature T(> T m ) 
and I < l c , the system can reduce the its free energy by 
decreasing the thickness of liquid film until the thickness 
reaches to zero, that is, the system crystallizes. When 
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I > l c , the system can reduce its free energy by increasing 
the thickness of liquid film until the thickness reach to 
infinite, the GB melts. l c is the critical thickness at T. 
If at a temperature l c = 0, any small thickness can lead 
to the melting of GB. The temperature is named as the 
maximum superheated temperature T Sl expressed as 



T„,a 



2Lpf 



(6) 



l c is also be expressed as 



T — T 

1 CT I s m \ 

l c = £Ln{ ) 

1 J-m. 



(7) 



The critical nucleus is extremely large for T w T m 
and reduces rapidly with the increase of the superheat- 
ing degree. At T s , l c = (Fig.^c)) and the spontaneous 
melting happens. Between T m and T s the superheated 
states is meta-stable, although the melting doesn't oc- 
cur. In following several sections, by MD simulations of 
Aluminum GB melting and above model, we prove that 
GB can preserve crystalline even above T m until temper- 
ature reaches to the maximum superheated temperature 
T s . Our simulation shows the behavior of superheating 
GB, that is, for T m <T < T s , there exists a critical width 
l c of liquid film. When the width of the artificially added 
liquid is larger than l c the GB melts, or the liquid film 
will reduce and the effect of crystallization is dominating. 



III. THE MOLECULAR SIMULATION OF 
GAIN BOUNDARY 
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FIG. 2: The distribution functions P(S) (a) and the pair cor- 
relation functions (b) at various temperatures. The flat and 
wide distribution of S(K) at 950K represents the coexistence 
of liquid and solid. 



Molecular Dynamic simulations of crystal 32] and 
high-symmetry surface |16| have shown the superheat- 
ing crystal without critical nuclei such as point defects 
and liquid drops. By internal heating, crystal with high- 
symmetry fcc(Ill) surface will be superheating. This is 
because the process of nucleation is homogeneous in the 
surface region. For some low symmetry surfaces such as 
fcc(llO) surface, the anisotropy leads to heterogeneous 
nucleation and high concentration of defects, superheat- 
ing is extremely difficult to achieve. 

Reconstruction and roughening arc two main struc- 
tural transitions for a surface when increasing temper- 
ature. Behavior of grain boundary is more complicated 
than surface, including the migration, bending, sliding, 
zigzag and faceting transition. Our simulations show that 
it is difficult to control the homogeneous nucleation in 
computer simulation to obtain superheated grain bound- 
ary. However we can use symmetric grain boundary to 
make atoms homogenously distribute in GB region. In- 
stead of periodic boundary, we have fixed two boundaries 
of simulation cell parallel to GB plane to decrease the 
possibility of grain-boundary sliding. 

In this work we choose Al symmetric £13 (320) [001] 
tilt grain boundary Fig^d) as our simulation cell with 



mis-orientation 67.8°, the tilt axis of the boundary is 
along the Z direction, the fixed boundary-conditions is 
used in the X direction perpendicular to GB plane, and 
periodic boundary condition are used in the Y and Z di- 
rections. The widths along X, Y, Z are 100A, 43l, 122 
respectively. Interatomic potential plays the most im- 
portant role in molecular dynamics simulations. In this 
paper, by using a more realistic potential, a Gl ue p oten- 
tial developed by F. Ercolessi and J. B. Adams [23], The 
glue potential has been used in a large number of the sim- 
ulations of surface, cluster, liquid, and crystal, and the 
simulation results are perfectly consistent with experi- 
mental results HE |3g. The lattice constant a 
for Al at OK is 4.032A. For the glue potential the melting 
points is about 936K (3^, which is close to the experi- 
mental melting point of Aluminum (about 933K). The 
MD simulations are carried out at constant temperature, 
constant volume and constant atomic number. The time 
step is 0.06 ( about 0.003psec ), at each temperature the 
runs are made about 20000 time steps ( about 60psec ). 

The static structure-factor 5(K) for specific reciprocal- 
space vector K (Eq.|8J represents a quantitative measure 
for the long-range order and can be used as order param- 
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(b) 1050K 
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(c) 975K (d ) 975K 

FIG. 3: The configuration of GB at high temperature, (a) The GB structure at 850K (high temperature disorder state) (b) 
The GB structure at 1050K (melting state) (c,d) The GB structure at 975K (superheated state) with unstable liquid film. 



eter describing the transition between disorder and order. 

S(K) =<| £ ex P ( lK ■ *(*)) ! 2> I N gb (8) 
iegb 

ri(t) is the position of ith atom at time t. < ■■• > 
is indicative of the time average. Nqb is the number 
of atoms in grain-boundary region. For a crystal the 
S'(K) is approaching 1, for liquid it approaching to 0. 
In order to study the stability and the nature of dis- 
ordered grain-boundary, we define a distribution func- 
tion of P(S) which is the statistics of the value of S(K) 
(K = |^(0,0,1)) of all molecular-dynamics time steps. 
Fig. |21 shows the distribution functions at several tem- 
peratures from 400K to 1050K. The centers of the peaks 
represent the degree of disorder and order, the widths of 
the peaks as a measure of the fluctuation of S(K) rep- 
resent the instability of the GB structure. From Fig. 
the peak is very narrow and the center of the peak is 
very near 1 in low temperature regime (T < 400K) , and 
this implies that the structure of grain boundary is very 
order and stable in low temperature regime; At 850K, 
the peak is board ( the minimum is at about 0.5 and the 
maximum is at about 0.9 ) and the center of the peak is 
much less than 1 , which implies that the GB structure 
is disordered in this temperature. At about 950K, the 
width of peak is extremely board (the minimum reaches 



0. 1 and the maximum still retains at about 0.9), the fact 
implies the structure of grain boundary is rather unsta- 
ble, sometime the structure of grain boundary is rather 
disorder like liquid because of the S(K) approaching 0, 
and sometimes it just likes a crystal structure with long 
range order because of the S'(K) probably approaching 

1. This shows that at this temperature the coexistence 
of liquid phase and solid phase is reached. Above re- 
sults show that the melting-point of this system is very 
close to 950K. In our simulation, the large S(K) fluctu- 
ation near 950K just indicates the signal of solid-liquid 
phase transition. At 1050K the peak moves to the left 
and become narrow. The position of the peak is close 
to zero, the GB is melting. Our results also show that 
the maximum superheated temperature T s is close to 
1050K. We also calculated the pair correlation functions 
■9( r ) = TirpN < Si^j fi( r ~~ r ij) > & t various temperatures 
which show the liquid behavior at 1050K. Fig.|H{a,b) il- 
lustrate the GB structures at 850K and 1050K. 

Our model shows that between T m and T s , GB will en- 
ter a new superheated state of GB. The new state is char- 
acterized by (1) the coexistence of liquid and solid; (2) 
the smaller size of liquid nucleus than that of the critical 
nucleus at that temperature prevents from the melting 
of GB although T > T m . Fig. Efc,d) show the compe- 
tition of liquid and solid phase at 975K in superheated 
state. Sometimes there exists a liquid-like layer in the 
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(a) 




S9& 






FIG. 4: The response of superheated state for adding liquid film. (a,b) at 850K, 21 =3uA. (c,d) at 975K, 21 =2uA. (e,f) at 
975K, 21 =3uA. 



grain-boundary region (c) but it is meta-stable and may 
disappear and crystalline phase re-emerges at following 
time steps (d). We will show that the superheated state 
(975K) is rather different from the high-temperature dis- 
ordered state at 850K. 

A liquid layer as the melting nucleus can be artificially 
added to the GB region by following methods: At a cer- 
tain temperature T, we sample some layers with width 21, 
the atoms in these layers are heated up to an appropri- 
ate temperature T a (> T) until a liquid layer forms. By 
allowing all atoms relaxation at temperature T again, 
we can obtain a new equilibrium structure at tempera- 
ture T. Fig. 21 shows that both the initial GB configu- 
rations Fig. Hfa,c,e) having already added a liquid film 
with widths 21 and the final equilibrium configurations 
Fig. H>,d,f) o at 975K and 850K respectively. At 850K, 
for 21 — 30A, the liquid layers disappears after the re- 



laxation about 50ps Fig. E{a,b). However at 975K for 
21 = 20A, the liquid layer disappears Fig. QJc,d)and for 
21 = 30A the layer of liquid is still existent Fig.0Je,f) af- 
ter relaxation about 50ps. Therefore, we can obtain 20A 
< 2l c < 30l at T=975K. 

Above results also show that the superheated state 
(975K) is very different in nature from the high tempera- 
ture disorder state (850K) . For high temperature disorder 
state, the liquid film can't induce the melting of GB, but 
for superheated GB state, the GB melts only when the 
width of liquid film is larger than a critical width l c . In 
order to define the correlation length £ and the thickness 
of liquid film, we calculate the atom density profile cor- 
responding Fig.@ff). Our results show that at 975K, the 
critical width 21 c is between 20 A and 30 A and £ = 10 A, 
thus 1.0 < Z c /£ < 1-5, which is consistent with the theory 
model Z c /£ = 1.4 with T m = 950K and T s = 105QK in 
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FIG. 5: The profile of potential and atomic density across 
the GB at 975 (K). 



Eq.0 



We must justify that pressure plays important roles 
for the superheated in experiments and computer simu- 
lations. In experiments, both coating with high-melting- 
point materials and heating internally induce the internal 
pressure in melting region. In our simulation also there 
is internal pressure in melting region because the size of 
simulation cell doesn't change companying with the in- 
creasing temperature. In this paper we only consider the 
superheating state due to the nucleation mechanism. The 
pressure mechanism and nucleation mechanism become 
intertwined and influence the melting of materials. Pres- 
sure lead to the increase of melting points. The melt- 
ing point is about 950K in our simulation and higher 
than the experimental melting points T exp —933K and 
Tgi ue =936K [3?1 | in simulation using glue potential with- 
out internal pressure. Pressure increases melting point 
less than T m -T giue ~950K-936K=14K. However the nu- 
cleation mechanism leads to the superheating about T s - 
T m -100K. 

The proximity effects of the superheated grains are 
also important to induce the superheated grain boundary 
sandwiched between two properly superheating grains. 
Because superheated state is meta-stable state in phase 
diagram, we don't expect it's long-live. The superheated 
materials melt by the nucleation mechanism or homoge- 
neously melt at higher temperature. 



IV. DISCUSSION AND CONCLUSION 

In summary, we study the possibility of superheated 
GB state by both a theoretical model and MD simulation. 
Our results indicate that we can obtain superheated grain 
boundary by having properly controlled homogeneously 
nucleating precession when increasing temperature. If 
there are liquid nuclei whose sizes are larger than a crit- 
ical size the superheated grain boundary melts. Or the 
grain boundary waits for homogeneously melting when 
temperature higher than maximum superheated temper- 
ature. 
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